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We consider the out-of-equilibrium evolution of a classical condensate field (j> = ($) and its 
quantum fluctuations for a $ 4 model in f + 1 dimensions with a double well potential. We use the 
two-particle point-irreducible (2PPI) formalism in the two-loop approximation. We compare our 
results to those obtained in the Hartree approximation, in the bare vertex approximation (BVA) 
and in the two-particle irreducible next-to-leading order large-N (2PI-1/N) approach, with thermal 
initial conditions. In the 2PPI scheme we find that the system tends to the symmetric configuration 
at late times, as expected in the absence of spontaneous symmetry breaking. 
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I. INTRODUCTION 

In recent years there was considerable activity in ex- 
tending the study of quantum field theory out of equi- 
librium beyond the leading order one-loop, large-TV and 
Hartree approximations H II II H 11 ill 

The simple $ 4 quantum field theory serves as a toy 
model for more realistic models with more degrees of free- 
dom. Such models are e.g. candidates for a description 
of the inflationary phase of the early universe. 

Recent numerical simulations 0,11,11 within differ- 
ent approximation schemes beyond leading order have 
lead to contradictory results for the phase structure of 
0(1) and O(N) $ 4 models in 1 + 1 dimensions. From 
finite temperature field theory it is known |To|. that in 
1 + 1 dimensions there should be no spontaneous symme- 
try breaking except at zero temperature; hence the ab- 
sence of spontaneous symmetry breaking is expected in 
nonequilibrium quantum field theory as well. While the 
Hartree approximation gives a first order phase transi- 
tion, the situation is different for various next-to-leading 
order approximation schemes. 

Cooper et al. have studied in detail the symmetry 
structure of the 0(1) model in 1 + 1 dimensions and the 
equilibration of the one- and two point correlation func- 
tions within the bare vertex approximation (BVA) and 
the two-particle irreducible next-to-leading order large- 
N (2PI-1/N) scheme. They found a second-order phase 
transition in the BVA and no spontaneous symmetry 
breaking in the 2PI-1/N approximation. Both approx- 
imation schemes are based on the so called Cornwall- 
Jackiw-Tomboulis (CJT) effective action ^lj- I n the 
CJT formalism the effective action is expanded in terms 
of two-particle irreducible (2PI) diagrams. The self en- 
ergy fulfills a Schwinger-Dyson equation and resums this 
2PI diagrams. A related approximation proposed by Ver- 
schelde and Coppens 0, is the two-particle point- 
irreducible (2PPI) formalism, where only local contribu- 
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tions to the self energy are considered. The one-loop 
2PPI approximation represents the Hartree approxima- 
tion. We have studied, in Ref. 8§, the 0(1) model within 
the two-loop 2PPI approximation and, with the param- 
eters and initial conditions chosen there, we found no 
spontaneous symmetry breaking. In order to have a di- 
rect comparison with the results of Cooper et al. we 
present here numerical simulations using their parame- 
ters and initial conditions. Besides the question of spon- 
taneous symmetry breaking it is interesting to compare 
the general features in the temporal evolution of all three 
approaches. 

The plan of this paper is as follows: In section ITTl we 
briefly summarize the model formulated in 8] and its 
extension to the case of an additional thermal initial state 
of quanta. In section IIIII we present the results of the 
numerical calculations. We end with some conclusions in 
section llVl 



II. THE MODEL 

The Lagrange density for the 4> 4 quantum field theory 
is defined by 
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(2.1) 



Here we will consider the case of a double well potential 
with /i 2 < 0. Within the 2PPI formalism the inverse 
Green's function G~ 1 (x,x') = i[D + M 2 (x)]S(x — x') re- 
mains local and is parametrized by a variational mass 
Ai 2 (t). The 2PPI effective action in terms of the varia- 
tional parameters <fi{t) and A4 2 (t) denotes 
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+ ±cf(x) + ±-(M 2 (x)-n 2 ) 2 

+r 2PPI [0,A4 2 ] . (2.2) 

The term r 2PPI [</!>, M 2 } contains all 2PPI vacuum dia- 
grams, i.e. all graphs that do not decay if two lines 
meeting at the same point are cut. 
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FIG. 1: We display the leading diagrams in the 2PPI action 
r 2PPI : (a) the bubble diagram; (b) the sunset diagram; solid 
lines: internal propagators; dashed lines: external fields (f>. 

By variation of the effective action r[^,A^ 2 ] with re- 
spect to M 2 one easily checks that M. 2 (x) satisfies the 
gap or Schwinger-Dyson equation 
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M 2 (x)=fi 2 + -^(x) + -A(x) , 
where the insertion A(x) has been defined as 
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This form of the effective action in terms of </> and M 2 
is analogous to the one with <p and A proposed by Ver- 
schelde and Coppens 0, . 

By truncation of the infinite series of 2PPf graphs in 
p2PPi one 0D -t; am g an approximation for the effective ac- 
tion. The two graphs depicted in Fig. ^ represent the 
relevant contributions to F 2PPI in the two-loop approx- 
imation, that will be used here. If only the one-loop 
bubble graph is included we refer to it as Hartree ap- 
proximation. 

In the following we will consider only spatially homo- 
geneous fields. Time integrations are understood to be 
carried out along a closed time path (CTP) [l^ IIH . 

The initial density matrix of the system is character- 
ized by a nonzero value 0(0) for the mean field and a 
distribution function n p for the quanta. We will as- 
sume a Bose-Einstein distributed initial state of quanta 
parametrized by 
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with uip = \J 'p 2 + m 2 ,, where m 2 , is defined below. The 
Green's function can therefore be written in terms of 
mode functions as 
1 
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(2n p + l)[f(t,p)t(t',p)e(t-t') 



-f(t',p)f*(t,p)Q(t'-t)] 



(2.6) 



The mode functions f(t,p) satisfy the differential equa- 
tion 



f(t,p)+[p 2 +M 2 (t)] /(t,p) = 



(2.7) 



The initial mass mo = Ai(0) has to be determined self 
consistently from the renormalized equation 
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with 8ui n = -^-lnAr. The mode functions at t = 

r nn 8-7T 771q 

satisfy f(0,p) = 1 and f(0,p) — —iu) v . 

The renormalized equations of motion in the two-loop 
approximation are 

= 4>(t) + M 2 m(t) - ^ 3 (t) + S(t) (2.9) 

M 2 {t) = v 2 + 6nL + ±^ 2 (t) + A^(t) + M 2 Ht)) 

(2.10) 

The finite one-loop part Ag^ has been defined as 

A&(*) = J ^[{2n p + l)\f{t,p)\ 2 -l} . (2.11) 



The two-loop contributions are given by 
S{t) = -i^ J dt'W) 

' 3 3 
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and 
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A^(t) = -A 2 / dt'<f>(t') / di"0(t") 
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x[G(t,t';p 3 )-G(t',t;p 3 )] (2.13) 
x [G(t',t";p 1 )G(t',t";p 2 )G(t,t"; P3 ) 
-G(t",t';pi)G(t",t';p 2 )G(t",t;p 3 )] . 

For the BVA and 2PI-1/N approximation Cooper et 
al. introduce the mass term x(t) as an auxiliary field. 
Naively, in our scheme this corresponds to the quantity 



x(t) = M 2 + ^ 2 n + ^[0 2 W + A< 1 n ) W 

= M 2 (t)^^ 2 \t) 



(2.14) 



where we have identified TrG with A^ . We note, how- 
ever, that the Green's functions have a different meaning 
in their and our schemes; this should be kept in consid- 
eration when comparing the numerical results. 



III. RESULTS 

The numerical implementation for the solution of the 
equations of motion [see Ea. (|2.9|l and (|2.10|) ] has been 
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described in detail in Ref. |8|. The equations of mo- 
tion are solved using a Runge-Kutta algorithm with a 
time discretization At = 0.001. The energy is conserved 
within five significant digits. We have chosen /i 2 = — 1, 
0(0) = 0.4, T = 0.1 and A equal to 3.0 and 21.3. This 
choice of parameters corresponds to the one of Cooper et 
al. 0. 

The time evolution of the mean field <f>{t") is presented 
in Fig. and Fig. [3^ and the time evolution of x(t) is 
displayed in Fig. 03 and Fig. These are compared 
with the results obtained in the Hartree approximation, 
and to those obtained by Cooper et al. for the BVA and 
the 2PI-1/N approximations. For the simulation with 
A = 21.3 we additionally present the time evolution of 
the effective mass Ai 2 (t) in Fig. 01 and the energy contri- 
butions in Fig. |SJ 
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2: Time evolution for the double well potential. Parame- 
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FIG. 

ters: fx 2 = 1, A = 3.0, 0(0) = 0.4, To = 0.1 (a) time evolution 
of the mean field <j>(t); (b) time evolution of the field xM> 
the solid lines relate to the two-loop 2PPI approximation, the 
dashed lines to the one-loop 2PPI or Hartree approximation, 
the dotted lines to the 2PI-1/N approximation and the dashed 
dotted lines to the BVA approximation, where the last two 
simulations are taken from |j| 

Surprisingly the mean field obtained in the 2PPI ap- 
proximation takes off rather early from the results ob- 
tained in the other approximations. It starts to oscillate 
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FIG. 3: Same as Fig. EJbut for A = 21.3 



around <f> = for t > 4, this indicates the absence of 
spontaneous symmetry breaking. The amplitude of the 
oscillations decreases with time. Remarkably, the mean 
field of the 2PI-1/N approximation decreases with the 
same time scale and similar amplitude, but without os- 
cillations, forming a kind of envelope. 

The composite field drives to negative values at 
early times and oscillates around positive values for later 
times [Figs. EJd and^]. For the simulation with A = 
3.0 the transition between the two regimes sets in at a 
time t > 13, whereas for the simulation with A = 21.3 
it takes place for t > 4. While in the simulation with 
A = 21.3 the average value of the field x(t) m the two- 
loop 2PPI approximation is, for later times, very close to 
the BVA and 2PI-1/N results, the situation for A = 3.0 is 
manifestly different. Here the average value of %(i) in the 
two-loop 2PPI approximation lies significantly below the 
one obtained in the other approximations. As we have 
mentioned in the previous section the comparison is not 
on firm grounds. 

Both the mean field <j> and the field x m the two- 
loop 2PPI approximation oscillate more strongly than 
the ones of the BVA and 2PI-1/N approximations. How- 
ever the amplitude is much smaller than the one of the 
Hartree approximation, and it decreases with time. 
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FIG. 4: Time evolution of the effective mass M 2 (t) for the 
simulation in Fig. the solid line represents the two-loop 
2PPI approximation, the dashed line the one-loop or Hartree 
approximation 
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FIG. 5: Time evolution of the energy for the simulation in 
Fig. [3 the solid line is the total energy, the dashed line the 
classical energy and the dotted line the quantum energy 



In Fig. 2] we display the time evolution of the vari- 
ational mass M 2 {t) in the two- loop 2PPI approxima- 
tion for the simulation with A = 21.3. As the differ- 
ence of x(t) and A4 2 (t) is directly proportional to A^ 2 ) (t) 
[see Eq. I2.14| , one concludes that a large contribution in 
Ai 2 (t) comes from the two-loop term [Fig. Ub]- The en- 
ergy displays strong oscillations at the transition towards 
the symmetric phase. At late times the energy has been 
transferred to the quanta almost entirely. 



IV. SUMMARY, CONCLUSIONS AND 
OUTLOOK 

We have compared numerical simulations in the two- 
loop 2PPI approximation to the Hartree approximation 
and to results in the BVA and 2PI-1/N approximation 
by Cooper et al.. 

In contrast to the Hartree and BVA approximations 
the two-loop 2PPI approximation exhibits no sponta- 
neous symmetry breaking in the noncquilibrium evolu- 
tion of $ quantum field theory. The mean field <p(t) 
tends to the symmetric configuration with = 0. The 
same is true for the 2PI-1/N approximation. It is some- 
what surprising that an approximation which is less pow- 
erful in rcsumming higher order graphs fares well in re- 
producing the expected phase structure. We think that 
this deserves further investigation. 

In 3 + 1 dimensions the two-loop 2PPI approximation 
has been applied 0,0] to the 0(1) and O(N) models in 
thermal equilibrium. There the approximation displays 
a second-order phase transition as expected for the exact 
theory. It therefore seems promising to study the out-of- 
equilibrium evolution in these models as well. Work on 
this is in progress. 

We thank Stefan Michalski for helpful conversations. 
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